Event-by-event study of neutron observables in spontaneous and thermal fission 
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The event-by-event fission model FREYA is extended to spontaneous fission of actinides and a 
variety of neutron observables are studied for spontaneous fission and fission induced by thermal 
neutrons with a view towards possible applications for detection of special nuclear materials. 



I. INTRODUCTION 

Phcnomenological studies of nuclear fission are of par- 
ticular interest for possible practical applications in the 
fields of nonproliferation and security. In particular, the 
detection of special nuclear material (SNM) has risen in 
priority. To better exploit all means of SNM detection, 
new efforts are underway to improve neutron detection 
technology, especially for the study of fast fission neu- 
trons from nuclear material. Since all SNM emits neu- 
trons, it is advantageous to use these neutron emissions 
for the detection of such material. For example, in highly 
enriched samples of plutonium (90% '^^^Pu, 10% 240pu) 



and uranium 



235 



U. 



10%238U), 



the small content of 



^^°Pu and 238u undergoes spontaneous fission, emitting 
on average two neutrons per fission. If it were possible to 
employ observable differences in the characteristics of the 
fission process between the two components of the mate- 
rial, it might be possible to distinguish between enriched 
and non-enriched samples of SNM. 

Being of a penetrating nature, neutrons may provide 
specific signatures of SNM and thus have intrinsic ben- 
efits over other observables. Their long attentuation 
lengths imply that they can propagate further through 
shielding material than photons emitted by fission. Un- 
fortunately their low cross sections in material constitute 
a natural drawback to neutron-based detection methods: 
the longer time required to obtain a clear measurement 
and the large solid angle necessary for detectors to sub- 
tend to collect as many neutrons as possible ^j]. 

The neutron background is relatively low, especially 
for prompt fission neutrons, emitted by the fragments di- 
rectly after fission. The dominant ambient neutron back- 
ground to fission neutrons is from cosmic rays. Neutron 
emission is thus a fairly unique signature of fissile mate- 
rial. Unfortunately, measurements of the neutron energy 
spectra alone are not particularly useful for determining 
the isotopic content of a sample material. Interactions in 
matter, particularly through a shielding material, shifts 
the scattered neutrons to lower energy, causing the yield 
of energetic fission neutrons to be reduced by several 
orders of magntiude relative to the peak of the distri- 
bution. Nevertheless, it is advantageous to study these 
neutrons using new fast-response detection systems be- 
cause fast neutrons preserve their direction better than 
thermal neutrons and can thus be used in SNM detection 



schemes involving correlated observables [l[. The present 
paper addresses such signatures of fission. 

Heretofore, most fission simulations have assumed that 
all emitted neutrons are drawn from the same energy 
spectrum which precludes correlations between the neu- 
tron mutliplicity and the associated spectral shape. In 
our event-by-event treatment, such inherent correlations 
are automatically included and we examine them with an 
eye toward specific applications. Our approach employs 
the fission model FREYA (Fission Reaction Event Yield 
Algorithm) which incorporates the relevant physics with 
a few key parameters determined by comparison to data 
• It simulates the entire fission process and produces 
complete fission events with full kinematic information on 
the emerging fission products and the emitted neutrons 
and photons, incorporating sequential neutron evapora- 
tion from the fission fragments. (We will examine prompt 
fission photon production in a later publication.) FREYA 
provides a means of using readily-measured observables 
to improve our understanding of the fission process and 
it is, therefore, a potentially powerful tool for bridging 
the gap between current microscopic models and impor- 
tant fission observables and for improving estimates of 
the fission characteristics important for applications. 

We compare and contrast correlations between neutron 
observables in spontaneous and thermal fission of ^'"'Pu, 
^'"^Pu(sf) and ^■^^Pu(nth,f) respectively, as well as be- 
tween ^3^U(sf) and ^3^U(nth,f). We also study these ob- 
servables in the spontaneous fission of ^^^Cf, often used as 
a calibrator for other fission measurements, and '^'^'^Cm. 

In the next section, we describe the experimental data 
employed, in particular the fission fragment mass distri- 
butions and the total fragment kinetic energy as a func- 
tion of fragment mass. We then discuss various neutron 
observables, including the prompt fission neutron multi- 
plicity as function of fragment mass, the neutron mul- 
tiplicity distribution, and the energy spectrum of the 
prompt fission neutrons. We also study the neutron- 
neutron angular correlations as well as the correlations 
between both the total kinetic energy of the fission prod- 
ucts and their residual excitation energy as a function of 
the total neutron multiplicity. We finally discuss poten- 
tial ways to exploit these correlations and conclude with 
some final remarks. 
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II. MASS AND CHARGE PARTITION 

The treatment of spontaneous fission in FREYA is simi- 
lar to that of neutron-induced fission, except for the sim- 
pUfication that there is no possibihty for any pre-fission 
emission. Thus, generahy, we start with a fissile nucleus 
"^"Zq having a specified excitation energy Eq, and let 
it undergo binary fission into a heavy fragment ^"Zh 
and a complementary light fragment ^^Z^. The frag- 
ment masses are obtained from experimental mass yields 
by the procedure employed in the original description of 
FREYA Q. 

The fragment mass yields, Y{A), are assumed to ex- 
hibit three distinct modes of Gaussian form [7] , 



Y{A) ^ SM) + S2{A) + Sl{A) 



(1) 



The first two terms represent asymmetric fission modes 
associated with the spherical shell closure at iV = 82 
and the deformed shell closure at iV = 88, respectively, 
while the last term represents a broad symmetric mode. 
The symmetric mode is relatively insignificant for spon- 
taneous fission which is at rather low nuclear excitation. 
The exception is ^^^Cf with a comparatively large sym- 
metric contribution. 

The asymmetric modes have a two-Gaussian form. 



N, 



27r(T,- 



(2) 



while the symmetric mode is given by a single Gaussian 



Sl - 



Nl 



,-{A-AYl2ai 



with A — ^Aq. Since each event leads to two fragments, 
the yields are normalized so that 'Y1,a ^(^) = 2. Thus, 



2iVi -f 27V2 + Nl 



(4) 



apart from a negligible correction because A is discrete 
and bounded from both below and above. 

The results are shown for the fission fragment and the 
subsequent product yields in Figs. [T][3] The yields re- 
ported for spontaneous fission are for the (primary) frag- 
ments while those reported for indued fission are for the 
(post-evaporation) products [S] . The modeling of the fis- 
sion fragment yields for neutron-induced fission over a 
range of incident neutron energies is discussed in Ref . Q . 
The product yields are obtained after FREYA has finished 
emitting neutrons from the excited fragments. All the 
yields exhibit similar behavior, a rather broad double- 
humped distribution with a gap near symmetry, A^jl. 
The symmetric contribution is typically very small. 

The results in Fig. [T] are most closely related because 
both ^■^^Pu(nth,f) and ^'*°Pu(sf) start from a compound 
nucleus with the same value of ^o- The ^'*°Pu(sf) data 
were taken from a study of 238,24o,242p^|^g£^ relative to 
^^^Pu(nth,f) The experiment was set up next to a 

reactor so that ^^^Pu(nth,f) could be used as a calibrator. 
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FIG. 1: (Color online) The percent yield as a function of 
fragment mass for 2='8pu(nth,f) ;8] (top) and ^''°Pu(sf) 
(bottom). The product yield data are shown for ^'^'■'Pu(nth,f) 
while the fragment yields are given for ^''''Pu(sf). The black 
curves are the 5-Gaussian fits to the fragment distributions 
while the red curves are the results after neutron emission in 
FREYA. 



with a large acceptance geometry to partially compen- 
sate for the low rate of spontaneous fission. (The high- 
est collected total number of spontaneous fission events, 
about 12000 for 240pu(sf), was a factor of about 200 be- 
low the number of thermal neutron-induced events.) The 
somewhat larger widths of the mass distributions result- 
ing from thermal neutron-induced fission were attributed 
to increased intrinsic excitation energy near the scission 
point 

The ^^^Cf(sf) fragment yields, shown in the upper 
panel of Fig.[2l result from an analysis of 2.5 x 10^ events 
[13| . The experiment focused on the far asymmetric mass 
region and showed that the enhancements in the yields 
observed previously were due to the choice of angular se- 
lection criteria. Choosing cos0 > 0.9 eliminated events 
where energy loss in the foil is large. We choose these 
results for use in FREYA, even though our focus is not 
on the far- asymmetric region, because the large sample 
size provides more accurate input. The ^^•^Cf(sf) yields 
were also measured with ^^*Cm(sf) in Ref. [Ijl, albeit 
with fewer statistics, 5.4 x 10^ ^^^Cf fissions and 71000 
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FIG. 2: (Color online) The percent yield as a function of 
fragment mass for 252Cf(sf) 1^ (top) and 2**Cm(sf) 14] (bot- 
tom). The data are fission fragment measurements. The black 
curves are the 5-Gaussian fits to the fragment distributions 
while the red curves are the results after neutron emission in 
FREYA. 



Cm fissions. The Cm fragment mass distributions 
and average fragment masses reported in Ref. are 
in good agreement with previous measurements, thus we 
can accept their reported yields with some confidence. 

The •^•^^U(sf) yield data in Fig. |3]was obtained from a 
uranium sample with a natural isotopic composition, i.e. 
with a small admixture of ^35^ 15]. A ^^^Cf(sf) neutron 
source was placed outside a double ionization chamber 
to provide thermal neutrons and thus allow comparison 
of the fission characteristics of ^^^U(sf) with ^^^U(nth,f). 
The double ionization chamber allowed measurements of 
the yields and kinetic energies of the two fission fragments 
in coincidence. They corrected for energy loss in the 
backing material and excluded angles greater than 60° 
to eliminate events where the fragments passed through 
more material, necessitating a larger correction for en- 
ergy loss in matter. The results shown here are based 
on 2800 fission events. They noted more fine structure 
in ^'^*U(sf) than in thermal neutron-induced fission of 
235^^ The yield at ^ ^ 119 is very poorly determined 
p^ . Unfortunately, no other ^•^*U(sf) data were found 
for comparison. 




A 

FIG. 3: (Color online) The percent yield as a function of frag- 
ment mass for ^3^U(nth,f) ;8;] (top) and ^^**U(sf) "15] (bottom). 
The product yield data are shown for ^'^^U(nth,f) while the 
fragment yields are given for ■^^*U(sf). The black curves are 
the 5-Gaussian fits to the fragment distributions while the red 
curves are the results after neutron emission in FREYA. 



In all cases, the locations of the asymmetric peaks in 
the data are similar while there appears to be a clearer 
separation of the asymmetric peaks at symmetry for the 
spontaneously fissioning nuclei. The asymmetric Gaus- 
sians also appear to be somewhat narrower in the case of 
spontaneous fission with the exception of ^'^^Cf, as shown 
in Fig. [2] In this case, the tails of the asymmetric dis- 
tributions shown are quite broad so the dip at symme- 
try is filled in to a considerable degree, even relative to 
neutron-induced fission. 

A clear shift between the fragment yields (before neu- 
tron emission) and the product yields (after neutron 
emission) is apparent in all cases. The magnitude of the 
shift depends on the overall mean neutron multiplicity, z7, 
which in turn depends on the partition of the excitation 
energy between the light and heavy fragments. The shift 
is not symmetric but is larger for the light fragment, es- 
pecially near symmetry. The location of the peak in the 
heavy fragment yield at A ^ 130 does not exhibit a sig- 
nificant shift due to neutron emission in any of the cases 
shown, even though the shift is apparent for other values 
of A. This is due to the proximity of the doubly-magic 
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closed shell with Zu — 50 and iV// — 82. This behavior 
is also apparent in the shape of TKE(A^f ) and in the de- 
pendence of the mean neutron multiplicity on fragment 
mass, as will be discussed later. 

The fragment charge, is selected subsequently. For 
this we follow Ref. Q and employ a Gaussian form. 



cx e 



(5) 



with the condition that \Zf — Zf{Af)\ < haz- The cen- 
troid is determined by requiring that the fragments have, 
on average, the same charge-to-mass ratio as the fission- 
ing nucleus, Zf{Af) = AfZo/Ao. The dispersion is the 
measured value, az = 0.5 @. The charge of the comple- 
mentary fragment then follows using Zl + Zh = Zq. 



III. FRAGMENT ENERGIES 

Once the partition of the total mass and charge among 
the two fragments has been selected, the Q value asso- 
ciated with that particular fission channel follows as the 
difference between the total mass of the fissioning nucleus 
and the ground-state masses of the two fragments, 



Qlh = M{A„)-Ml-Mh 



(6) 



FREYA takes the required nuclear ground-state masses 
from the compilation by Audi et al. [l^, supplemented 
by the calculated masses of Moller et al. [11| when no 
data are available. The Qlh value for the selected fis- 
sion channel is then divided up between the total kinetic 
energy (TKE) and the total excitation energy (TXE) of 
the two fragments. 

Figures |4][6] show the measured average TKE value as 
a function of the mass number of the heavy fragment. 
Ah . Near symmetry, the plutonium fission fragments are 
mid-shell nuclei subject to strong deformations. Thus 
the scission configuration will contain significant defor- 
mation energy and the TKE will be correspondingly low. 
At Ah — 132, the heavy fragment is close to the doubly- 
magic closed shell having Zh — 50 and Nh — 82 and 
is therefore resistant to distortions away from spheric- 
ity. Consequently, the scission configuration is fairly 
compact, causing the TKE to exhibit a maximum even 
though the complementary light fragment is far from a 
closed shell and hence significantly deformed. Note that 
the peak around Ah = 132 is a feature of all the data 
sets shown, regardless of whether fission is neutron in- 
duced or spontaneous and independent of the identity of 
the fissile nucleus. 

The ^'^^Pu(nth,f) data sets in the top panel of Fig. |3] 
are very consistent for Ah > 135, above the closed shell 
at Ah = 132. In this region and below, the agreement 
among the data sets is not as good, particularly near 
the symmetry value of Ah = 120, presumably due to 
the low fragment yields in this region. Unfortunately, no 
uncertainties are given on the data, only the full-width 
half maximum spread of TKE for several given values of 
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FIG. 4: (Color online) The total fragment kinetic energies as 
a function of the heavy fragment mass for 239pu(nth,f) [13- 
[Hl (top) and ^*''Pu(sf) [12] (bottom). The FREYA results are 
shown with the calculated variance arising from the range of 
charges available for each Ah- 



Ah in the measurment of Nishio et al. [l7|. This vari- 
ance is similar to that shown for FREYA. The ^'"'Pu(sf) 
data by Schillebeeckx et al. [ijj are somewhat flatter in 
the region of the closed shell. Unfortunately there are 
considerable fluctuations in the data for Ah < 130 and 
TKE(yl//)was not measured for Ah < 122. We have 
therefore extrapolated a constant average value back to 
Ah = 120. 

The ^^^Cf(sf) data in Fig. [5] are again taken from 
Ref. [11] with cos 6* > 0.9. The high statistics of this 
measurement result in small experimental uncertainties 
and smooth behavior of TKFi{Ah). There is more un- 
certainty in the lower statistics ^''^Cm(sf) data from 
Ref. hj. Indeed, a comparison of earlier measurements 
in that work showed that although the average fragment 
masses were consistent, the average TKE of ^'*^Cm(sf) 
varied by 4% among measurements, depending on the 
measurement techniques as well as the choice of calibra- 
tors (either ^^^Cf(sf) or ^^^U(nt,h,f)). The results shown 
in the bottom half of Fig. [5] agree with the highest re- 
ported energy and indeed are ~ 2 MeV higher than those 
reported from ^^^Cf(sf) [l3]- No more recent results on 



5 



210 
200 




150 - 
210 — 



2001- 
^ 190h 



^ ISOh-r 
X 

< 

m 170 
160 
150 



252 



Cf(sf) 



244 



* FREYA ^ 
o Schmitt-Henschel 



Cm(sf) 



J 



120 



130 



140 150 

Ah 



160 



170 



190 



^ 180 

> 





170 



^ 160 



> 

CD 



150 
190 

180 

170 



< 

^ 160 



150 



235 



U(n,f) 



o Nishio 



», 



238 



U(sf) 



* FREYA 
o Ivanov 



120 130 140 150 160 



H 



FIG. 5: (Color online) The total fragment kinetic energies as 
a function of the heavy fragment mass for ■^^■^Cf(sf) [31 (top) 
and ^*'*Cm(sf) [31 (bottom). The FREYA results are shown 
with the calculated variance arising from the range of charges 
available for each Ah- 



FIG. 6: (Color online) The total fragment kinetic energies 
as a function of the heavy fragment mass for ■^'^^U(nth,f) [20] 
(top) and ^^*U(sf) ^ (bottom). The FREYA results are shown 
with the calculated variance arising from the range of charges 
available for each Ah- 



^'*''Cm(sf) are available. 

There are significant fluctuations in the ^'^^U(sf) data 
[Tsj in Fig. |6| which can be attributed to the rather low 
statistics of this measurement. We note that the data in 
Ref. [l^ were presented as a function of the light frag- 
ment mass instead of the heavy. The distribution shown 
here is obtained by reflection. We note also that the 
symmetry region is rather poorly measured with points 
missing around Ah ~ 119. The remaining data are at 
values of TKE below the lower limit of the plot. Per- 
haps some of the fluctuations in TKE can be attributed 
to the 'fine structure' noted in the yields in Ref. ^15] . 
They also noted that their measured average TKE was 
on average 3 MeV lower than their calibrated result for 
thermal neutron-induced fission of ^'^'^U. 

We assume the average TKE values take the form 

TKE(Aff,i;n) = TKEdata(AH)+dTKE(£;„) . (7) 

The first term on the right-hand side of Eq. ([7|) is ex- 
tracted from the data shown in Figs. HKH while the sec- 
ond term is a parameter adjusted to ensure reproduction 
of the measured average neutron multiplicity, V. In each 
particular event, the actual TKE value is then obtained 



by adding a thermal fluctuation to the above average, as 
explained later. 

Figures llHni include the average TKE values calculated 
with FREYA at thermal energies for neutron-induced fis- 
sion and for spontaneous fission, together with the as- 
sociated dispersions. Thus the bars associated with the 
FREYA calculations are not sampling errors but indicate 
the actual width of the TKE distribution for each Ah- 

IV. NEUTRON EMISSION 

Once the average total fragment kinetic energy has 
been obtained, the average combined excitation energy 
in the two fragments follows by energy conservation, 

TXE = El + Eh ^ Qlh- TKE - (8) 

The first relation indicates that the total excitation en- 
ergy is partitioned between the two fragments. As is 
common, we assume that the fragment level densities are 
of the form Pi{E*) ~ exp{2 \/aiUi), where Ui is the effec- 
tive statistical energy in the fragment and Oi is the level- 
density parameter. We follow the prescription of Ref. [3| 
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FIG. 7: (Color online) The distribution of the maximum tem- 
perature in the residual nucleus, P{T), after v neutrons have 
been evaporated from a primary fragment for ^■'''Pu(nth,f) 
(top) and 2''°Pu(sf) (bottom). The curves show results for 
v = 1 (dashed), u = 2 (dot-dashed), and f = 3 (dot-dot- 
dashed) as well as the distribution including all i/ (solid). 



FIG. 8: (Color online) The distribution of the maximum tem- 
perature in the residual nucleus, P{T), after v neutrons have 
been evaporated from a primary fragment for ^^^Cf(sf) (top) 
and ^**Cm(sf) (bottom). The curves show results for = 1 
(dashed), v = 2 (dot-dashed), and u = 3 (dot-dot-dashed) as 
well as the distribution including all v (solid). 



with the value of the asymptotic level density parameter 
Co obtained from the ^■^^Pu evaluation, assuming it to be 
universal. 

If the two fragments are in mutual thermal equilib- 
ritim, Tl =Th, the total excitation energy will be pro- 
portional to the level-density parameters, i.e. ^ ai. 
FREYA therefore first assigns tentative average excitations 
based on such an equipartition, 

E* = ^^li^L^TXE, (9) 

where E* = {Ai/Ao)TXE. Subsequently, because the ob- 
served neutron multiplicities suggest that the light frag- 
ments tends to be disproportionately excited, the average 
values are adjusted in favor of the light fragment, 

eI^xEI ,Eh = tke-eI , (10) 

where x is an adjustable model parameter expected be 
larger than unity. 

After the mean excitation energies have been assigned, 
FREYA considers the effect of thermal fluctuations. The 



fragment temperature Ti is obtained from Ui = Ui{E*) — 
aiTf and the associated variance in the excitation E* is 
taken as = 2U*Ti, where U{E*) = E* in the simple 
(unshifted) scenario. 

Therefore, for each of the two fragments, we sample a 
thermal energy fluctuation SE* from a normal distribu- 
tion of variance af and modify the fragment excitations 
accordingly, so that 

E* - E:+SE: ,z = L,H. (11) 

Due to energy conservation, there is a compensating op- 
posite fluctuation in the total kinetic energy Q. The 
corresponding dispersions are included in Figs. |4]|6l 



A. Temperature distributions 

FREYA assumes that the (fully accelerated) fission frag- 
ments first deexcite by (possibly sequential) neutron 
evaporation, followed by sequential photon emission. At 
each stage in the decay chain, the spectral shape of the 
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FIG. 9: (Color online) The distribution of the maximum 
temperature in the residual nucleus, P{T), after v neutrons 
have been evaporated from a primary fragment for ^^"'U(nth,f) 
(top) and ^^**U(sf) (bottom). The curves show results for 
V = 1 (dashed), v = 2 (dot-dashed), and v = i (dot-dot- 
dashed) as well as the distribution including all (solid). 



ejectile is determined by the maximum temperature in 
the daughter nucleus, which in turn is determined by 
the excitation in the emitting nucleus and the associated 
Q- value. Because there are many different Q- values in- 
volved, one for each fragment species, and becaue the 
excitation in the emitter fluctuates, the maximum tem- 
perature in the daughter nucleus displays a non-trivial 
distribution. 

The distributions of the maximum temperature in the 
evaporation daughter, granddaughter, and great grand- 
daughter nuclei are shown in Figs. [7][n] for the various 
cases under study. The maximum temperature in the 
daughter nuclei (^"=1) is fairly well peaked at around 
0.5-0.8 MeV; the larger contribution comes from the light 
fragment which tends to be hotter than its heavy partner. 
The temperature distributions in the granddaughter nu- 
clei {v—'i) are considerably broader and peak at lower en- 
ergies, In the case of the great-granddaughters (t^=3) the 
distributions decrease monotonically. For spontaneous 
fission of ^"^^Pu and ^^^U, where the average total neu- 
tron multiplicity is only V is 2.15 and 2.0, respectively, 
it is most likely that each fragment emits just a single 




Fragment mass A 

FIG. 10: (Color online) The neutron multiplicity as a func- 
tion of fragment mass for neutron-induced fission of ^^®Pu 
[itI, [isl, l2l|] (top) and spontaneous fission of ^''"Pu (bottom). 
The vertical bars at intervals of A = 10 on the data from 
Ref. [l7j in the top panel are the full-width half maximum 
of the distribution in neutron multiplicity from the range of 
charges available for that A. The FREYA results also show this 
variance by the vertical bars. (Note that the scales on the 
y-axes of the plots are not identical.) 



neutron so the probability that a fragment emits three 
neutrons is rather small, and the probability for further 
neutron emission is, in most cases, too small to be visi- 
ble on the plot, though such events do contribute to the 
overall distribution near zero. 

The induced fission processes ^■^^Pu(nth,f) and 
^■^'^U(nth,f) lead to higher daughter temperatures than 
the spontaneous processes ^^''Pu(sf) and ^'^*U(sf) and 
they also yield larger multiplicities. Thus 77 « 2.88 for 
239pu(nth,f) relative to 2.15 for 240pu(sf), and v « 2.47 
for 235u(nth,f) relative to 2.0 for 238u(sf). Relative to 
^'"^Pu(sf) and ^■^*U(sf), spontaneous fission of the heavier 
actinides, ^^"'Cm and ^^^Cf, also results in higher daugh- 
ter temperatures. Indeed P{T) is similar for ^'*^Cm(sf) 
and ^'^^Pu(nth,f), as are the average neutron multiplic- 
ities, 2.72 and 2.88, respectively. We note that for 
^^^Cf(sf) the maximum daughter temperature peaks at 
T 1 MeV, which is higher than in all the other cases 
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Fragment mass A 

FIG. 11: (Color online) The neutron multiplicity as a function 
of fragment mass for spontaneous fission of ^^^Cf 23- 25] (top) 
and ^'''*Cm 22] (bottom). In the bottom part of the plot, the 
points labeled uncorrected are also from Ref. [2^ but do not 
include corrections for fragment motion. The FREYA results 
show the variance in ^{A) from the range of charges in the 
vertical bars. (Note that the scales on the j/-axes of the plots 
are not identical.) 



considered, and the tail extends up to T « 2 MeV. 
In addition, the distribution of the maximum temper- 
ature in the granddaughter has a distinct peak around 
T « 0.6 MeV, which is not surprising considering that 
i7w 3.75 for 252cf(sf). 

In the treatment by Madland and Nix [l^ a convenient 
analytical expression was obtained by assuming that the 
overall distribution of the maximum daughter temper- 
ature has a triangular shape. The results in Figs. [7][9] 
demonstrate that such an assumption is hardly realistic, 
even for just v = 1. 



B. Average neutron multiplicity 

The dependence of the average neutron multiplicity on 
the fragment mass number A, is very sensitive to the di- 
vision of the excitation energy which is governed by the 
parameter x in Eq. (ITUl) . As shown in Figs. [imiT^ all 




Fragment mass A 

FIG. 12: (Color online) The neutron multiplicity as a function 
of fragment mass for neutron-induced fission of ^^^U 20, 26] 
(top) and spontaneous fission of ^'^*U. The FREYA results show 
the variance in ^{A) from the range of charges in the vertical 
bars. (Note that the scales on the j/-axes of the plots are not 
identical.) 



the measurements exhibit a characteristic 'sawtooth' be- 
havior: the neutron multiplicity from the light fragment 
increases slowly as A approaches \Aq and then drops 
rather sharply to a minimum around Ah ^ 130, the same 
location as the maximum of TKE(Ah) ■ Due to the pres- 
ence of the closed shell at that point, the fragments are 
particularly resistent to neutron emission. Past the dip 
region, the multiplicity again increases. The dip tends 
to be more sharply defined for larger nuclei where ^Aq is 
close to 130. For example, the drop is particularly abrupt 
for ^^^Cf where ^Aq — 126. Where data are available, it 
is seen that the FREYA calculations provide a rather good 
representation of the 'sawtooth' behavior of 'S'(^), even 
though FREYA is not tuned to these data. 

Although the agreement is generally good, the ob- 
served behavior is not perfectly reproduced. The FREYA 
results for ^^^Pu(nth,f) in Fig. [TOl agree very well with the 
data for 90 < A < 140 with a; = 1.1. At higher and lower 
A, although there are deviations, the measurements are 
within the variance of the FREYA results. However, in 
these regions, as well as near symmetry, the yields are 
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Neutron multiplicity v 

FIG. 13: (Color online) The probability for a given neu- 
tron multiplicity as a function of multiplicity, v, for neutron- 
induced fission of ^'^^Pu [i^l (top) and spontaneous fission 
of ^""Pu 28, 29] (bottom). The FREYA results are shown by 
the black circles while the equivalent Poisoon distribution is 
shown by the red squares. 



smaller so that larger deviations may be expected. No 
neutron measurements were made in Ref. |12| and we 
have not located any comparison data for ^'*''Pu(sf); here 
we use X = 1.2. Below the symmetry point, the slope of 
T'{A) is rather small, but it increases more rapidly above 
A — 132. The FREYA results display larger fluctuations 
for this isotope, which may be attributed to the more ir- 
regular behavior of TKE(A^f). Also note that while the 
central values of ^^{A) fluctuate more, the variances are 
the same size as for ^■^^Pu(nth,f) in the top panel and 
appear exaggerated by the smaller scale of ^{A) on the 
y-axis for ^'"'Pu(sf). 

There are numerous measurements of ^{A) for 
^^^Cf(sf). A sample of some representative, more recent, 
results are shown in the top part of Fig. 1111 The mea- 
surements are all very similar with small differences only 
near A < 90, A 120, and A > 150. The light fragment 
data are rather flat and then increase rather quickly for 
105 < A < 120. Above A = 132, the slope oil7{A) is less 
than that seen for 239Pu(nth,f) in Fig. \M The FREYA 
results, obtained with x — 1.3, while consistent within 




Neutron multiplicity 



FIG. 14: (Color online) The probability for a given neutron 
multiplicity as a function of multiplicity, y, for spontaneous 
fission of ^'^^Cf m (top) and ^'^''Cm (bottom). The 

FREYA results are shown by the black circles while the equiv- 
alent Poisoon distribution is shown by the red squares. 



the variance of the data, are significantly flatter than the 
data for A < 100 and A > 140. This behavior, stronger 
than any of the other FREYA results in Figs. fimiT^ can be 
traced to the apparent two-slope behavior of TKF,{Ah) 
for ^^^Cf(sf) in Fig. [5l there is a slow decrease in TKE 
for 132 < Ah < 145 with a faster change of TKE with 
Ah thereafter. It is also consistent with the large widths 
of the asymmetric fission yields in Fig. [51 We note that 
changing x does not change the slopes of i^iA) , only the 
relative magnitudes. To better describe T'{A) with FREYA, 
it would be necessary to fit x{A) rather than employing 
just a single- valued parameter for x. 

Data for ^'^*Cm(sf) are shown in the lower panel of 
Fig. [Til B oth data sets shown are from the same experi- 
ment [2 2( 1 . Two results were given in the paper due to the 
location of the fission source relative to the neutron de- 
tector, the corrected version, labeled Schmidt-Henschel, 
shifts the measured F(A) to account for the fact that only 
some of the neutrons emitted by the fragments will reach 
the detector. The authors calibrated their correction for 
^^^Cf(sf) by normalizing the position of the sawtooth 
with A to previously published data. The ^^'*Cm(sf) cor- 
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rection was made by scaling the ^^^Cf shift by the ratio of 
the total neutron multiplicity in ^'*'*Cm relative to ^^^Cf, 
presumably because they could not compare their ^'^^Gm 
results to other measurements of the same system. The 
authors were rather inconclusive about which results were 
actually correct since the sum of complementary multi- 
plicities, 17(A) + ^{Ao — A), did not agree well with the 
total neutron multiplcity, Ft (A) near ^^o- If this dis- 
agreement is real, then the correction had the effect of 
shifting the peak of 17 (A) downward from A '-^ 120 to 
A ^ 117. Interestingly, the FREYA results, calculated with 
X — 1.2 and treating the motion of both the fragments 
and the neutrons relativisticall y, ag rees much better with 
the uncorrected results of Ref. [22l ■ 

Finally, FREYA results are compared to ^^^U(nth,f) in 
the top panel of Fig. [TH The values of Vj A) agree very 
well with the sawtooth pattern of the data [20, [2^ with 
the exception of the symmetric region where the yields 
are rather low. Indeed, Ref. [20,] does not provide re- 
sults for 17(A) in the symmetric region. The FREYA re- 
sults for ^•^®U(sf), calculated with x — 1.2, are shown in 
the bottom panel of the figure. The shape of the saw- 
tooth appears rather flat for this isotope, likely because 
of the high 17(^A) obtained near symmetry, reflecting the 
low TKE reported by Ref. [15|] in Fig. H 



C. Neutron multiplicity distribution 




Neutron multiplicity v 



Figures 11311151 show the neutron multiplicity distribu- 
tion P(i^) for the various cases considered. Each emitted 
neutron reduces the excitation energy in the residue by 
not only its kinetic energy (recall E = 2T where T is the 
maximum temperature in the daughter nucleus) but also 
by the separation energy Sn (which is generally signifi- 
cantly larger). Therefore the resulting P(zy) is narrower 
than a Poisson distribution with the same average mul- 
tiplicity, as clearly seen in the figures. 

In experiments, the quantity Pii') is determined by 
detecting fission events in a sample of material and cor- 
relating these with simultaneous neutron detection. The 
relative probability for emission of v neutrons in given 
event, P{v), is inferred by combining the calculated prob- 
ability for observing n neutrons when u were emitted, 
Q{n] v), with the detector efficiency determined from the 
count rate by comparison with a calibration source hav- 
ing a known 77; typically ^^^Cf(sf) is used. Thus, while 
the value of V may be well measured for a given isotope, 
the distribution P{i') is less well determined. 

We compare to data in so far as possible. The results 
labeled 'Holden-Zucker' in Figs. [H] for 239pu(nth,f) and 
[15] for ^'^^U(sf) are consensus values from a 1985 report 
by Holden and Zucker [27j. Results from this reference 
are generally available for the other cases shown here. 
However, we do not show them if there is good agree- 
ment between the data displayed here and in Ref. [27| or 
if there are more recent data available. In most cases, 
the agreement is rather good, with the FREYA results fol- 



FIG. 15: (Color online) The probability for a given neu- 
tron multiplicity as a function of multiplicity, u, for neutron- 
induced fission of ^"^^U [s^-fs^ (top) and spontaneous fission of 
23*U 27] (bottom). The FREYA results are shown by the black 
circles while the equivalent Poisoon distribution is shown by 
the red squares. 



lowing the data more closely than the equivalent Poisson 
distribution. 

The largest difference between the calculated FREYA 
multiplicity distributions and the data seems to be for 
^'"^Pu(sf), which may be due to the smaller sample of 
240pu(sf) Y{A) and TKE{Ah) data used as input to 
FREYA. By contrast, the agreement of FREYA with the 
^■^^U(sf) consensus data given in Ref. [27'] is remarkable. 



D. Multiplicity-gated neutron spectra 

For fisison events having a specified total neutron mul- 
tiplicity v, we define the associated spectral shape, 



fniE) 



1 diy 



(12) 



which is thus normalized to unity, while the correspond- 
ing spectral shape of the neutrons from all the fission 
events irrespective of the associated multiplicity is de- 
noted simply by fniE) and is also normalized to unity. 
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FIG. 16: (Color online) The spectral shapes of prompt fis- 
sion neutrons for neutron-induced fission of ^^®Pu (top) and 
spontaneous fission of ^''"Pu (bottom), for events with a fixed 
neutron multiplicity of = 1, . . . , 6 as well as for all events 
irrespective of the associated multiplicity. 



FIG. 17: (Color online) The spectral shapes of prompt fission 
neutrons for spontaneous fission of ^^^Cf (top) and ^''''Cm 
(bottom), for events with a fixed neutron multiplicity of — 
1, . . . , 6 as well as for all events irrespective of the associated 
multiplicity. 



The multiplicity-gated spectral shapes obtained for the 
various cases considered are shown in Figs. [TMTBl Re- 
sults are presented for multiplicities up to = 6. It is 
apparent that the spectra become progressively softer at 
higher multiplicities, as one would expect because more 
neutrons are sharing the available energy. This type of 
elementary conservation-based correlation feature is not 
provided by the standard models of fission. 

The tails of the prompt fission neutron spectra from 
^'*''Pu(sf) are longer and broader than those from 
^'^^Pu(nth,f) even though the average energies are smaller 
and fewer neutrons are emitted. The opposite is the 
case for ^'^*U(sf) and ^''^U(nth,f) since the spectra from 
^^®U(sf) are closely clustered around the mean. The most 
energetic neutrons at high multiplicity are emitted from 
^^^Cf(sf) where the spectra are also rather closely clus- 
tered around the mean. The spectral shapes and average 
energies of ^^'*Cm(sf)are rather similar to ^'^^Pu(nth,f)- 

Table U shows the mean kinetic energy of neutrons 
emitted from the two fragments as a function of the 
neutron multiplicity. The averages are shown for multi- 
plicities with event samples large enough to prevent the 
results from being skewed by fluctuations. The largest 



average energies are seen for ^^^Cf(sf) with similar re- 
sults for ^^^Pu(nth,f) and ^'^'*Cm(sf), not surprising given 
the temperature profiles seen in Figs. [7]|ni The variances 
are on the same order as the averages, albeit somewhat 
smaller, in most cases. 



E. Neutron-neutron angular correlations 

The event-by-event nature of FREYA makes it straight- 
forward to extract the angular correlation between two 
evaporated neutrons, an observable that has long been of 
experimental interest (see, for example, Rcfs. [33, 36, 3'^ 
and references therein) but which cannot be addressed 
with the standard models of fission. 

Figures [TUll^ show this quantity for the neutrons re- 
sulting from fission induced by thermal neutrons on ^ssy 
and ^"^^Pu as well as neutron correlations in spontaneous 
fission. The results are shown for neutrons with kinetic 
energies above thresholds a.t E — 0.5, 1 and 1.5 MeV. The 
angular modulation grows somewhat more pronounced 
as the threshold is raised (while the statistics are corre- 
spondingly reduced). 
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FIG. 18: (Color online) The spectral shapes of prompt fis- 
sion neutrons for neutron-induced fission of ^^s-q (top) and 
spontaneous fission of ^-^^U (bottom), for events with a fixed 
neutron multiplicity of = 1, . . . , 6 as well as for all events 
irrespective of the associated multiplicity. 
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TABLE I: The mean neutron kinetic energy, (E), together 
with the associated dispersion, as, for events with a fixed 
neutron multiplicity v as well as for all events. 



The neutrons tend to be either forward or backward 
correlated. The backward correlation appears to be 
somewhat favored. While not shown, we have previously 
analyzed the case of ^■^^Pu(nth,f) for — 2, breaking 
it down to three separate contributions: both neutrons 
from the Hght fragment, both from the heavy fragment, 
and one neutron emitted from each fragment [sl] . There 
is a significant correlation at 9i2 = when both neutrons 
are emitted from the same fragment, with a higher peak 
for the case when both neutrons are emitted from the 
light fragment due to its higher velocity. On the other 
hand, when one neutron is emitted from each fragment, 
their direction tends to be anti-correlated due to the rela- 
tive motion of the emitting fragments, resulting in a peak 
at 9i2 ~ 180. The overall result is a stronger backward 
correlation because emission from both fragments is most 
likely. 

The backward correlation is strongest when the overall 
neutron multiplicity is low, especially for ^^"Pu(sf) and 
^'^*U(sf), whereas large multiplicities, as for ^^^Cf(sf) and 
^'*''Cm(sf), reduce the angular correlation. 



F. Correlations between product energies and 
neutron multiplicity 

The combined kinetic energy of the two resulting (post- 
evaporation) product nuclei is shown as a function of the 
neutron multiplicity ly in the top panels of Figs. 22-24. 
It decreases with increasing multiplicity, as one might 
expect on the grounds that the emission of more neu- 
trons tends to carry off more initial excitation energy, 
thus leaving less available for the products. As expected 
from the behavior of Z^Zh, the combined product ki- 
netic energy is largest for the most massive fission sys- 
tems (^^^Cf and ^''^Cm) and lowest for the least massive 
(235U and 238U). 

The bottom panels of Figs. 22-24 show the neutron 
multiplicity dependence of the average residual excitation 
energy in those post-evaporation product nuclei. Because 
this energy is available for subsequent photon emission, 
one may expect that the resulting photon multiplicity 
would display a qualitatively similar behavior and thus, 
in particular, be anti-correlated with the neutron multi- 
plicity. 
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FIG. 19: (Color online) The angular correlation between two 
neutrons emitted from ^^^Pu(nn,f) (top) and ^''°Pu(sf) (bot- 
tom) as a function of the opening angle between the two neu- 
trons, 6i2- The FREYA results are shown for several cuts on the 
neutron kinetic energy: E > 0.5 MeV (solid black), 1 MeV 
(dashed red), and 1.5 MeV (dot-dashed green). 



FIG. 20: (Color online) The angular correlation between two 
neutrons emitted in spontaneous fission of ^^^Cf (top) and 
^^^Cm (bottom) as a function of the opening angle between 
the two neutrons, The FREYA results are shown for sev- 
eral cuts on the neutron kinetic energy: E > 0.5 MeV (solid 
black), 1 MeV (dashed red), and 1.5 MeV (dot-dashed green). 



There is little sensitivity of the residual excitation to 
the identity of the fissioning nucleus in any of the cases 
presented. This result shows that the energies left over 
after prompt neutron emission are not strongly depen- 
dent on the initial fragment temperature. 

V. APPLICATIONS 

We have so far shown that there are characteristic cor- 
relations between the emitted neutrons that depend on 
relative angle, energy and multiplicity. To best take ad- 
vantage of these correlations, fast response detector sys- 
tems are desirable. Such systems can better exploit these 
correlations which would be washed out in slow response 
detectors and detection systems based on moderators. 
Scintillator detectors can distinguish between neutrons 
and photons with good background rejection. They are 
also amenable to scaling to larger solid angle coverage. 
Threshold detectors that gate on higher-energy prompt 
neutron emission can be readily used to study angular- 
energy or multiplicity-energy correlations [l[. 



As we have shown in Figs. [TM?T1 the neutron-neutron 
angular correlations can distinguish between configura- 
tions where both neutrons are emitted from a single 
fragment or one neutron is emitted from each fragment. 
This correlation will become stronger with neutron en- 
ergy, particularly for ^'^°Pu(sf) and ^^^U(sf) where the 
mean neutron multiplicity is rather low and the emit- 
ted spectra are softer. Such evident directionality could 
improve background rejection of neutrons from cosmo- 
genic sources. Comparison of correlations in admixtures 
of plutonium or uranium isotopes could reveal the degree 
to which the material is enriched. 

Figures fTBllTSl clearlv show the difference in the spectral 
shapes for specified neutron multiplicities. The slope of 
the spectral shape increases with neutron multiplicity so, 
for example, there are fewer energetic neutrons for 1/ = 3 
than for 1/ — 1. Again the difference in the spectral 
distributions gated on neutron multiplicity in neutron- 
induced fission relative to spontaneous fission in the same 
sample of material could be exploited by fast detector 
systems, providing an additional means of determining 
the isotopic content of the material. 
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FIG. 21: (Color online) The angular correlation between two 
neutrons emitted from ^^^U(nn,f) (top) and ^^*U(sf) (bot- 
tom) as a function of the opening angle between the two neu- 
trons, 012. The FREYA results are shown for several cuts on the 
neutron kinetic energy: E > 0.5 MeV (solid black), 1 MeV 
(dashed red), and 1.5 MeV (dot-dashed green). 



For experimental groups to better explore the possible 
correlation studies available with FREYA, we are providing 
a version to work in-line with several larger Monte Carlo 
codes, including MCNP [39]. More details about the in-line 
version of FREYA will be discussed in the future. 



VI. CONCLUDING REMARKS 

We have shown that event-by-event models of fission, 
such as FREYA, provide a powerful tool for studying fission 
neutron correlations. Our results demonstrate that these 
correlations are significant and exhibit a dependence on 
the fissioning nucleus. 

Since our method is phenomenological in nature, good 
input data are especially important. Some of the mea- 
surements employed in FREYA are rather old and statis- 
tics limited. It would be useful to repeat some of these 
studies with modern detector techniques. In addition, 
most experiments made to date have not made simul- 



taneous measurements of the fission products and the 
prompt observables, such as neutrons and photons. Such 
data, while obviously more challenging to obtain, would 
be valauble for achieving a more complete understanding 
of the fission process. 
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FIG. 22: (Color online) The total product kinetic energy (top) 
and residual excitation energy (bottom) remaining after neu- 
tron emission has ceased as a function of neutron multiplicity. 
The FREYA results are shown for ^^®Pu(nth,f) (squares) and 
2^°Pu(sf) (diamonds). 
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FIG. 23: (Color online) The total product kinetic energy 
(top) and residual excitation energy (bottom) remaining 
after neutron emission has ceased as a function of neutron 
multiplicity. The FREYA results are shown for spontaneous 
fission of ^''^Cm (squares) and ^^^Cf (diamonds). 



FIG. 24: (Color online) The total product kinetic energy 
(top) and residual excitation energy (bottom) remaining 
after neutron emission has ceased as a function of neutron 
multiplicity. The FREYA results are shown for ^^^U(nth,f) 
(squares) and ^■^**U(sf) (diamonds). 
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